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ABSTRACT 

Supernova remnants (SNRs) have long been assumed to be the source of cosmic rays (CRs) up to 
the "knee" of the CR spectrum at 10^^ eV, accelerating particles to relativistic energies in their blast 
waves by the process of diffusive shock acceleration (DSA). Since cosmic ray nuclei do not radiate 
efficiently, their presence must be inferred indirectly. Previous theoretical calculations and X-ray 
observations show that CR acceleration modifies significantly the structure of the SNR and greatly 
amplifies the interstellar magnetic field. We present new, deep X-ray observations of the remnant 
of Tycho's supernova (SN 1572, henceforth Tycho), which reveal a previously unknown, strikingly 
ordered pattern of non-thermal high-emissivity stripes in the projected interior of the remnant, with 
spacing that corresponds to the gyroradii of 10^^ — 10^^ eV protons. Spectroscopy of the stripes shows 
the plasma to be highly turbulent on the (smaller) scale of the Larmor radii of TeV energy electrons. 
Models of the shock amplification of magnetic fields produce structure on the scale of the gyroradius 
of the highest energy CRs present, but they do not predict the highly-ordered pattern we observe. We 
interpret the stripes as evidence for acceleration of particles to near the knee of the CR spectrum in 
regions of enhanced magnetic turbulence, while the observed highly ordered pattern of these features 
provides a new challenge to models of DSA. 

Subject headings: acceleration of particles — Cosmic rays — ISM: individual objects (SN 1572, Tycho) 
— ISM: supernova remnants 



1. INTRODUCTION 

Energetic arguments have long-favored the hypothe- 
sis that supernova shocks are responsible for the bulk 
of the cosmic ray (CR) flux at Earth from 10^ eV to 
the "knee" , a flatteni ng of the CR spec trum at approx- 
imately 10^^ eV (Bla ndford fc Eichlerj[l987| . This re- 
quires that ~ 10% of the mechanical energy of Galactic 
supernova remnants (SNRs) be converted into relativis- 
tic particles over the remnants' lifetime. In the favored 
model of this process, diffusive shock acceleration (DSA), 
self-generated magnetic turbulence provides scattering 
centers from which suprathermal particles may diffuse 
between the up- and downstream fluid , gai ning energy 
(on average) with each shock crossing ( [Bell |1978[ ). The 
ubiquitous non-thermal radio emission in SnKs is ev- 
idence of GeV electrons in SNRs, and the X-ray syn- 
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chrotron emission observed from the limbs of young SNRs 



CKoyama et al.'1995! 'H ^hes et al.|2QQQa[ [Gotthelf et al 
2001; Hwang et al. 2002| indicates the presence of elec- 
trons to tens of TeVs. Though CR nuclei comprise the 
vast majority of mass in SNR relativistic particles, they 
do not radiate efficiently, and their presence must be in- 
ferred indi rectly, either by the modification of the shock 
dynamic s (Warren et al. 2005 ) and thermal state of the 
plasma (Hughes et al. "2000b Decourchelle et al. 2000 
Patnaude et al. 20091), o r by tne produ ction of secondary 
p ions ( Butt et al. 2001; Aharonian et al. 2006). 
Bell & Lucek (2001 J predicted significant amplification 
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of the magnetic held in supernova shocks by CR- induced 
turbulent processes, allowing for acceleration to the knee 
within a SNR lifetime. This picture is supported by high 
spatial resolution observations from the Chandra X-ray 
Observatory that require high magnetic fields to explain 
the short synchrotron cooling time imp hed by the geo- 
metrically thin rim s of young remnants (Vink & Laming 
2003 Ballet 2006). An extensive theoretical literature 
has developed detailing a range of plasma instabilities 
that might generate these high fields. Of these, the cos- 
mic ray current driven instability (C RCD, ofte n referred 
to as Bell's non-resonant mechanism, |Bell|2004 ) , in which 
the magnetic field is amplified by t ur bulence generated 
by CRs propagating ahead of the shock, has received the 
most study, though other processes (e.g., acoustic insta- 
bilities, resonant instabilities, etc.) have been suggested. 
MHD tBeh 2005^ and particle -in-ceh (PIC) simulations 
|Kiquelme & ^pitkovskyj[2009|) show that the CRCD in- 



st ability drives plasma motions which suppress growth on 
small scales, until the dominant scale of the turbulence 
is comparable to the Larmor radius of the highest energy 
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Fig. 1. — Chandra X-ray 4.0-6.0 keV image of the Tycho supernova remnant, smoothed with a ~ 0.75^^ Gaussian and displayed with 
an arcsinh scahng, showing various regions of striping in the nonthermal emission. Clockwise from the upper right: a) The main western 
stripes discussed in this Letter; b) A fainter ensemble of stripes; c) a previously-known bright arc of non-thermal emission, with our newly 
discovered streamers; d) filaments of "rippled sheet" morphology common in optical observations of middle-aged SNRs. 



cosmic rays. This results in a no n- uniform medium with 
structured, filamentary density and field enhancements 
~ 10 times greater than the mean quantities, surround- 
ing cavities of size approximately equal to the Larmor 
radius of the highest energy particles present. Indeed, 
Bell's analysis and simulations show that once the insta- 
bility has saturated, the medium consists of a network of 
holes evacuated of plasma and field, enclosed by "wan- 
dering filaments" of high density and frozen-in field. One 
goal of observations is to discover evidence for structure 
in SNR shocks on this spatial scale. 

2. OBSERVATIONS 

We observed Tycho in April 2009 with the Chandra 
X-ray Observatory Advanced CCD Imaging Spectrome- 
ter imaging array, as part of a Cycle 10 Large Program 
(LP), using the 4 ACIS-I front-side illuminated CCDs, 
operated in Faint mode. The total program was split into 
nine individual obsids, which we reprocessed with CIAO 
version 4.1, using the gain tables and CTI correction in 
CALDB 4.1.3. Examination of the lightcurves revealed 
no significant background fiares. The registration of the 
individual pointings was improved using the measured 
relative positions of background point sources, though 
the necessary shifts were typically <Cl'^ For imaging 
analyses, the registered event lists were merged using the 
standard CIAO tools, and have a total average livetime 
of ~ 734.1 ks. For spectroscopy, counts were extracted 
and RMFs and ARFs were generated for each individual 
obsid, which we fit jointly in XSPEC version 12.6.0. 

In Figure [T] we show the image from the 4-6 keV band, 
which is dominated by the synchrotron component of 
the X-ray spectrum. Apart from the well-known limb- 



brightened shell, a number of bright regions are seen to- 
ward the projected interior of the remnant. While the 
strongest of th ese feature s are visible in earlier Chandra 
data (War ren et al.|2005 ) our deeper observation reveals 
a striking pattern or nearly-regularly spaced stripes. The 
brightest group, centered ~ 30'^ interior to the west- 
ern limb (see Figure [2]), has a peak surface brightness 
twice that of the brightest sections of the rims, and is 
the primary subject of this communication. A second, 
fainter pattern extends east- west, 55 — 75'^ inside the 
southern rim, and there is evidence for several other re- 
gions of striping near the detection limit. Comparison 
with a shallower 2004 Chandra observation reveals no 
statistically significant change in the brightness of the 
stripes, ruling out any dramatic fiux variability like that 
observed in the no n-thermal X-ray filame nts of the SNR 
RX J1713. 7-3946 ( |Uchiyama et al.||2007D . Our prelim- 
inary proper motion measurements tor the stripes are 
consistent with the overall expansion of the blast wave 
and, in particular, show no evidence for non-radial fiow. 
There are no obvious counterp arts to these features in 
the radio ( Reynoso et al.||1997 ), nor in the mid-IR. 



3. ANALYSIS 

3.1. Location of the Stripes 

Before investigating the nature of the stripes, we first 
must locate them within the three-dimensional volume 
of the remnant. The canonical picture of a young SNR 
consists of three distinct fiuid discontinuities: the blast 
wave, which marks the shock propagating into the am- 
bient medium, a Ray leigh- Taylor (R-T) unstable con- 
tact discontinuity (CD) at the ejecta-interstellar material 
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Fig. 2. — (Left) The projected intensity profile for two regions containing the western stripes. The bright stripe is the red peak around a 
position of 35^^ The locations of prominent stripes are indicated by the horizontal lines with tick marks. The peaks in the northern profile 
are spaced by 9.9'', 11.3'', 6.9", 8.3", and 13.3"; those in the southern profile are spaced by 4.4", 9.4", 6.8", and 7.4". (Right) Grayscale 
image of the Chandra data in a line- free region of the spectrum (4.2-6.0 keV) showing the extraction regions for the profiles plotted in the 
left panel. The profiles run from the lower left to the upper right and the tick marks on the boxes correspond to the same tick marks as 
in the left panel. The scale bar in the lower right is 10" long. The green boxes denote the spectral regions analysed in the Letter. The 
brightest region is the bright stripe; the coaddition of the fainter regions constitutes the faint stripes. 



boundary, and a reverse shock that p ropagates into the 
stehar remains. Warren et al. (2005) set an upper-hmit 
for the azimuthahy-dependent projected radius of the re- 
verse shock in Tycho using the location of the Fe Ka 
emission. Adopting their center of expansion, the west- 
ern stripes peak at a radius of 220'^ weh outside the 190'^ 
position of the reverse shock at that azimuth. While the 
position of the stripes does coincide with the Warren et 
al. estimate of the CD, the regularly- spaced, linear mor- 
phology of the non-thermal stripes does not correspond 
to any features in the R-T plumes of thermal emission 
tracing the ejecta boundary, nor is the CD a prominent 
feature elsewhere in the 4-6 keV band. Conversely, the 
blast wave is a bright source of 4-6 keV emission, and we 
identify the stripes as projected features of this forward 
shock. 

Tycho 's blast wave is traced by a very thin shell of X- 
ray em ission, with a typica l thickness only 1-2% of its 
radius (Warren et al. 2005). Since the stripes are seen 
in projection away from the rim, their line-of-sight path- 
length through the shell is small. Thus their intrinsic 
emissivity must be high relative to the limb regions. We 
examined two regions of different brightness on the west- 
ern rim with minimal ejecta contamination, and assumed 
a spherical shell model with a radius of 214^', chosen to 
match the curvature of the western limb. The radial 
brightness profile of the brighter rim region is too nar- 
row and peaked to be consistent with a simple projected 
shell geometry, meaning that its high surface brightness 
likely arises from a local enhancement in emissivity, per- 
haps similar to the stripes. In contrast, the profile of the 
fainter region is well- fit by a shell with thickness 2% of 
its radius, and is likely more typical of the blast wave 
as a whole. This simple model predicts that features at 
80% of the shock radius (where the stripes are located) 
would have a surface brightness 20% that of the limb. 



The observed peak surface brightness of the stripes is, 
however, 5 times that of the faint limb region. Thus 
the peak intrinsic emissivity of the synchrotron-emitting 
plasma in the stripes must be a factor of 25 higher than 
is typical for the blast wave. In principle, this brighten- 
ing could be due to a local enhancement in the ambient 
density and/or magnetic field. While we can not en- 
tirely exclude this possibility, we consider the existence 
of a pre-existing structure around Tycho with the right 
combination of increased density, magnetic field, and, as 
we show below, turbulence necessary to produce the ob- 
served correlations between the ordered structure and 
spectral variations improbable. In contrast, models of 
CR-driven magnetic field amplification produce structure 
in the precursor density and field, and we find it more 
likely that the stripes mark a region extensively modified 
by the acceleration process. This assumption is implicit 
in the analysis in this Letter. 

3.2. Spectroscopy 

The spatial position of the stripes, seen in projection 
against the bright thermal X-rays of the shocked ejecta, 
greatly complicates spectroscopy. This can be mitigated 
by restricting analysis to the 4.2-10 keV range, where 
the emission is dominated by the synchrotron continuum 
and the Fe Ka line complex at ~6.45 keV. Throughout 
this work we set the absorbing column to Njj = 7 x 
10^^ cm~^, which has negligible effect at these energies. 
The spectrum of the brightest stripe is well-fit with a 
power-law of photon index of F = 2.1lto.io5 significantly 
harder than the spectrum of the entire SNR at these 
energies (F = 2.72 ± 0.02; this fit includes faint lines 
identified by Tamagawa et al. 2009|), and a coaddition 



of several of the lower surface brightness stripes (F = 
2.52 ± 0.05; this fit includes an Fe K line). 
While photon indices are useful for characterization 
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Fig. 3. — (Left) Broadband spectra of Tycho from the radio to the X-ray for the entire SNR (top), an ensemble of faint stripes (middle) 
and the bright stripe (bottom). The radio points for the latter two came from our own estimates based on the VLA radio data, taking 
account of the remnant's expansion. The integrated radio flux density values for the entire SNR are taken from Klein et al. (1979). The 
Chandra data are plotted in red and for the top two spectra the Fe K line is visible. The overplotted curves are the best-flt phenomenological 
model for the continuum emission of a power law (extrapolated from the radio band) with an exponential cut-off in the X-ray regime. 
(Right) Chandra spectra of the three spatial regions analyzed. The top part shows the observed spectral intensity versus photon energy 
(data points with error bars) along with the best flt model. The bottom part shows the residuals (data — model) in units of the flux error. 
The data points and model curves here include instrumental effects, which have been removed for the X-ray data points plotted on the left 
panel. 



of the non-thermal spectrum, the broadband spectrum 
graduahy departs from a simple power law over several 
decades in energy. Detailed DSA models of synchrotron 
spectra show a small degree of curvature through the ra- 
dio to UV, with a cut-off in the X-rays, set (in young 
SNRs like Tycho) by the shor t cooling time of TeV-scale 
electrons (Ellison et al. 2000). Absent detailed spectral 
calculations tor 'lycho, we adopt a common approxima- 
tion and assume a power-law electron momentum dis- 
tribution with an exponential cut-off, which produces a 
photon spectrum that is nearly a power-law fr om the ra- 



dio through the UV and is curved in X-rays (Reynolds 



fc Keohane|[l999 ). For this model, the low energy index 
(ck^radio) and the flux normalization are set by radio ob- 
servations, while the location of the high energy cut-off is 
constrained with the X-ray spectrum. For the spectrum 
of the entire remnant, assuming a radio spectral index of 
Q^radio = 0. 65, a flux density at 1 GHz of ^50 Jy (Kothes 



et al [2006 ) , and again inclu ding the faint X-ray emission 
lines from Tamagawa et al.| ([2009j), the fit to the Chan- 
dra data gives z/cut = 1-9 x 10^ Hz {hucut ^ 0.79 keV). 
This provides a remarkably good fit {xl = 1.05), and is 
consistent with measurements from Suzaku. (The statis- 
tical error on this number is so small as to be irrelevant; 
the precision to which we measure it is far greater than 
the efficacy of our emission model or the accuracy of the 
radio flux and spectral index measurement s . ) If we in- 
stead assume aradio = 0-61 (|Klein et al.||1979|), the fit 
is slightly worse {xl 
lower energies {vcut — 



= 1.10), and the cut-off moves to 
0.9 X 10^^ Hz ; hucut 0.37 keV). 
For the stripes, we estimated the radio flux in these 
regions from a VLA map. The fits favor slightly shal- 
lower indices (ofradio ~ 0.60) and higher cut-off energies. 



z^cut = I9I9 X 10 Hz (7.91^;^ keV) for the Bright Stripe 
and z/cut = 2.8 ± 0.2 x 10^^ Hz (1.16 ± 0.04 keV) for the 
ensemble of fainter stripes. We plot the broadband syn- 
chrotron spectra of these regions in Figure |3] 

Although it is clear that the Bright Stripe and Faint 
Ensemble spectra require different best-fit parameter val- 
ues, in principle these differences might be due to either 



different cut-off frequencies or radio spectral indices. Our 
spectral analysis (detailed in Table |l]) allows us to dis- 
criminate between these two options. When the spec- 
tra for the Bright and Faint regions are fit with radio 
flux densities fixed at their nominal values, both regions 
prefer an index very close to the average for the entire 
remnant. Conversely, if we instead require that both re- 
gions have the same cut-off, the fitted radio indices are 
a = 0.565 (bright) and a = 0.623 (faint). We cannot 
reject this level of variation outright since measurements 
of the spatial variation of the radio spectral index are not 
available for Tycho, though in Cas A (a plausible com- 
parison object) the radio index d oes vary spatially by 
much as ±0.1 (^Wrig ht et al.||1999[). On the other hand, 
the best fit for a single cut-off frequency is clearly worse 
(Ax^ = 13.8) than for independent cut-off values. An 
F-test shows that the fit including independent cut-off 
frequencies is an improvement at a significance level of 
more than 3a. Thus we conclude that the cut-off fre- 
quencies are significantly different in the bright and faint 
spectral regions analysed. 

4. INTERPRETATION 

4.1. The Diffusion Coefficient 

DSA requires magnetic field inhomogeneities strong 
enough that particles can scatter multiple times across 
the shock front. This diffusion of CRs in the magnetized 
turbulence is generally assumed to be close to the Bohm 
regime, where the mean-free path is the particle gyrora- 
dius. Departure from this limit can be parameterized as 
ko = 1^0 /^o, Bohm, where I^o is the diffusion coefficient 
at the electron cut-off energy. Parizot et al. (2006) show 
that for an electron distribution limited by synchrotron 
losses, the photon cut-off energy {E^^cut = ^^cut) is set 
by the diffusion coefficient, shock velocity (V3, in units 
of 1000 km s~^), and compression ratio, r: 



0.14 



7,cut,keV *^ 3, shock 



16 (r- 1) 
372 



(1) 



At a distance of 4.0 ±1.0 kpc ( [Hayato et aL]|2010D, with 
a blast wave proper motion of 0.3'' yr~^ (Katsuda et al. 
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2010), Vs = 5700 kms ^. For an unmodified shock (r = 
4) and our fit cut-off energies, ko = 5.8 (Whole SNR), 
3.9±0.1 (faint stripes), and O.dto.t (Bright Stripe). With 
a higher compression more apphcable to Tycho (r = 7), 
these numbers become 3.8, 2.6 ± 0.1, and 0.4^^;^. (The 
error bars reflect only the statistical error in Vcut^ and 
not the ^ 25% uncertainty in the distance to Tycho.) 
Our value for the whole SNR is near to those found by 
[Parizot et al.,(2006) for Tycho (^o = 4.9 - 10), while the 
very small numbers for the brightest stripe are similar 
to their results for G347.3-0.5 and SN 1006 {ko = 0.2 - 
0.9). Since ko > I for isotropic turbulence, this may 
be due to our oversimplified emission model, and could 
indicate that the photon spectrum falls more gradually to 
high_energies, as seen in some DSA calculations (Ellison 
et al.j2010|), or may be indicative of anisotropic diffusion 



locity of 5700 km s ^, fits to models of Tych o 's syn- 



( [Keville et al. 2008). Regardless, our measured photon 



cut-ott energies require an order of magnitude decrease 
in the diffusion coefficient in the stripes relative to the 
average value in Tycho, and indicate that diffusion in 
the stripes is very near the Bohm limit. Furthermore, 
since Bohm diffusion requires 6B/B ^ 1, ko ^ 1 implies 
fully developed turbulence, at least on the scale of the 
gyroradius of the highest energy (~ TeV) electrons in 
the brightest stripes. 

4.2. The Origin of the Stripes 

Might the spatial structure and enhanced turbulence 
in the stripes be a manifestation of the plasma instability 
driving the cosmic ray acceleration? The largest possible 
characteristic scale of the acceleration process is set by 
the gyroradius of the highest energy particles present. If 
we identify the gaps between the stripes (/gap ~ 8^') with 
twice the proton gyroradius, an estimate of the magnetic 
field yields a measurement of the energy, given by Equa- 
tion [1 

The choice of B is somewhat uncertain. If the stripes 
are a consequence of the highest energy CRs interact- 
ing with the non-amplified ISM field {B ^ 3 /iG), for 
a gap spacing of 8'^ Ecr = 2 x 10^^ eV. However, 
these features may arise from a region where the field 
has already been somewhat amplified. For a shock ve- 



chrotron emission ( [Cassam-Chenai et al. 2007) predict 
an upstream field strength of ~ 30 /^G, indicating 
Ecr = 2 X 10^^ eV — just at the knee of the Galactic 
CR spectrum. 

5. CONCLUSION 

If the Bell mechanism is active in the magnetic field 
amplification in SNR shocks, and if the pattern of stripes 
we observe in Tycho is an observable consequence of this 
process, the stripes are direct evidence of particles ac- 
celerated to the knee. There are several caveats to con- 
sider before accepting this conclusion. Most obviously, 
while MHD and PIC simulations of the CRCD generate 
structure in the amplified medium, they do not show the 
regular and nearly periodic pattern that we observe in 
Tycho. Secondly, when the CR back-reaction is included 
(Riquelme fc Spitkovsk y|20"Q9 ), PIC simulations indicate 
that the Bell instability alone may be neither fast nor 
efficient enough to provide the observed amplification, 
particularly for Tycho's low ambient density. Additional 
instabilities are required, and seem to be indicated by 
the degree of turbulence in the stripes impl ied by our 



observations. Riquelme fc Spitkovsky] (2010) propose a 
new instability (the Parallel Current Driven Instability, 
PCDI) that acts on smaller scales (i.e., lower energies) 
and provides significant field amplification over that es- 
tablished by the CRCD. Further theoretical work on the 
degree of turbulence generated by proposed mechanisms 
on TeV particle scales seems particularly warranted. Re- 
gardless, the match in scales between the observed spac- 
ing of the Tycho stripes and the Larmor radius of knee- 
region CRs — a scale intrinsic to the CRCD — is tan- 
talizing, while the ordered structure presents a challenge 
to current models of diffusive shock acceleration. 
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TABLE 1 
Stripe Spectroscopic Fits 



Region 


a 


(1 GHz) 

(Jy) 


Z^cut 

(Hz) 


x2/d.0.f. 


Bright Stripe 


0.61* 


3.33 X 10-2* 


2.60 ±0.09 X lO^s 


194.89/205 


Bright Stripe 




3.33 X 10-2* 


1.88i J-j^ X 1018 


194.30/204 


Bright Stripe 




6.66 X 10-2* 


1.93+2-09 X 1018 


194.26/204 


Bright Stripe 


0.65* 


8.0+2-3 X 10-2 


2.37+i-g^ X 1018 


194.27/204 


Bright Stripe 


0.61* 


A.2ti-\ X 10-2 


1.561 j;^^ X 1018 


194.27/204 


Faint Ensemble 


0.61* 


0.39* 


3.52 ±0.03 X lOi'^ 


667.29/687 


Faint Ensemble 


n «ni +0.009 


0.39* 


2.84 ±0.10 X lOi'^ 


666.31/686 


Faint Ensemble 




0.78* 


2.88 ±0.11 X lOi^ 


666.55/686 


Faint Ensemble 


0.65* 


899+°-°^2 


3.53 ±0.23 X lOi^ 


666.37/686 


Faint Ensemble 
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0.394 ±0.012 


3.47i°:232 X 101^ 


667.25/686 
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